Photon-number entangled states generated in Kerr media with optical parametric 

pumping 



A. Kowalewska- Kudlaszylfl 

Nonlinear Optics Division, Department of Physics, 
Adam Mickiewicz University, Umultowska 85, 61-614 Poznan, Poland 

W. LeonskQ 

Quantum Optics and Engineering Division, Institute of Physics, 
University of Zielona Gora, Prof. Z. Szafrana 4a, 65-516 Zielona Gora, Poland 



Jan Pefina Jr. 
Palacky University, RCPTM, Joint Laboratory of Optics, 
17. listopadu 12, 771 46 Olomouc, Czech Republi^ 
^ ■ (Dated: January 20, 2013) 

, ' Two nonlinear Kerr oscillators mutually coupled by parametric pumping are studied as a source 

Oh| of states entangled in photon numbers. Temporal evolution of entanglement quantified by negativity 

. shows the effects of sudden death and birth of entanglement. Entanglement is preserved even in 

asymptotic states under certain conditions. The role of reservoirs at finite temperature in entan- 
glement evolution is elucidated. Relation between generation of entangled states and violation of 
Cauchy-Schwartz inequality for oscillator intensities is found. 



^ (— l PACS numbers: 03.67.Bg, 03.65.Yz, 42.50.Ex, 42.65.Lm 




* | ann akowQamu . edu. pi | 

t wlconski@proton. if . uz . zgora. pi 
t pcrinaj@prfnw.upol.cz 



2 



I. INTRODUCTION 

Generation of entangled states and their protection from disentanglement is a crucial task in quantum information 
processing theory. Entangled states can be generated in various physical systems. Entanglement can be observed in 
various degrees of freedom. The most common systems considered as sources of entangled states are spin systems 
(NMR), quantum dotsQ, trapped ions Q, double- well Bose-Einstein condensates 0). Typical constituents of these 
systems can be represented by two (or more) level atoms that mutually interact and create this way entangled states. 
We speak about systems composed of qubits (see for example [|| and the references quoted therein) , qutrits or higher 
dimensional qudits. On the other hand, entangled states can be found also in field systems, namely in optical fields. 
Entanglement emerges here due to mutual interaction among modes of optical fields. This interaction typically occurs 
in nonlinear optical processes [||. Entanglement can then be observed in polarization states (for example [gj]) or 
photon-number states 0- Namely states entangled in photon numbers are promising for the near future due to the 
progress in construction of photon-number resolving detectors in the last years. Their construction can be based 
on optical-fiber loop interferometers (8l-[l0|. absorption in superconductor wires 11 1, use of intensified CCD cameras 
[l~2l fl3| , complex semiconductor structures [l4[ or hybrid photodetectors dH • 

Every real physical system interacts with its environment [1(1] . Such interaction leads to losses of coherence, and 
consequently, to gradual degradation of entanglement. In quantum information processing this represents one of 
the crucial problems. That is why many methods for coping with this problem have been developed; quantum 
error correction [l7l - l20j decoherence free subspaces pl| - |23| , proper preparation of qubit states [24| , multiple system- 
decomposition method [25( , to name few. It is well known that an entangled state interacting with a zero temperature 
Markovian environment asymptotically decays in time. However the character of zero-temperature reservoir decay can 
be altered. Under specific conditions disentanglement may suddenly occur (sudden death of entanglement) 26 28]. 
Moreover, if a system interacts with a thermal or squeezed reservoir the entanglement is lost in finite time as 
well. We should also mention that the effect of sudden entanglement birth can also be observed [30M32| . 

In this paper, we consider two optical modes generated by parametric down-conversion [33| pumped by an intense 
optical field. The generated two optical modes propagate through a medium characterized by Kerr nonlinearity that 
modifies their statistical properties. States characterizing two optical modes can be entangled in photon numbers 
(see for example @) as a consequence of the fact that photons are emitted in pairs in parametric down-conversion. 
Provided that interaction with external reservoir is taken into account the effects of sudden death and sudden birth 
of entanglement can be observed. Entanglement will be quantified by negativity function and its presence will be 
compared with nonclassical behavior of distributions of integrated intensities of two modes. We note that experimental 
characterization of nonclassical two- mode optical fields is experimentally available [HJ . 

Optical fibers arc suitable candidates for experimental generation of the predicted entangled states. They allow 
both the generation of photon pairs in parametric processes (38l . [3^ and exhibit Kerr nonlinearities (37j . Nonlinear 
fibers allow to generate sufficiently high photon-pair fluxes. For example, a nonlinear fiber 300 m long embedded 
into a Sagnac interferometer allows to generate up to 10 7 photon pairs per 1 mW of optical pumping [38| . Similar 
photon-pair generation rates have also been reported from a 2 m long micro-structured fiber [39(. 

However, we note that nonlinear systems excited directly by an external coherent field [referred as nonlinear quantum 
scissors |40h42| |] require giant Kerr nonlinearities that can be achieved, e.g., using induced transparency [43j-|45|. 

The paper is organized as follows. A quantum model of two interacting optical fields is formulated and solved in 
Sec. II that also contains considerations about entanglement. Generalization of the model to the case of interaction 
with reservoir is provided in Sec. III. Nonclassical behavior of integrated intensities of two modes is discussed in 
Sec. IV. Sec. V provides conclusions. 

II. THE MODEL, ITS SOLUTION AND ENTANGLEMENT EVOLUTION 

The considered system is composed of nonlinear medium (characterized by Kerr nonlinearity) that is excited by an 
external field of strength g. At quantum level, two photons (belonging to different field modes) are emitted together 
in the parametric process. Interaction Hamiltonian Hi nt of the system can be written in the following form: 

Hint = #Kcrr + -ffpar (1) 

£ Kerr =f (aW + f (W (2 ) 

tfpar - gaW + g*ab. (3) 

Here, Hamiltonian H^ crr describes two nonlinear Kerr oscillators (with nonlinearities Xa,b) whereas Hamiltonian H pa[ 
means an effective interaction Hamiltonian of the two- mode parametric process [46| . Model Hamiltonians including 
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both parametric and Kerr processes have been analysed in [47|. It should be noted that our model does not include 
the Kerr cross term a'ab'b, that is important in some Kerr media (for discussion concerning models involving such 
Kerr cross term see, for example [48l |49| and the references quoted therein). We note that the Hamiltonian analogous 



to that considered here has also been used for example, by Bernstein |50j or Chcflcs and Barnett |51| . The dynamics 
of the models excluding the Kerr cross term have also been analyzed in other papers, for instance in [52l [53l|. 

It should be stressed out that various interesting effects beside of those discussed in the quantum information 
theory can originate in systems endowed with Kerr nonlinearities. For instance, Miranowicz at al. [54j have shown 
that discrete superpositions of an arbitrary number of coherent states (so-called Schrodinger cat-like or kitten states) 
can be generated in these systems. Moreover, such nonlinear systems (especially, nonlinear couplers with Kerr 
nonlinearities (40l - |4^ |) can serve as a source of maximally entan gled states (MES) useful not only in quantum optics 
but also in other branches of physics. For instance, Kurpas et al. [55[ discussed nonlinear resonances and entanglement 
generation in solid state models, particularly in those involving SQUID systems. 

At the beginning we neglect damping. That is why we simply apply the Schrodinger equation that can be written 
in the form of a set of equations for complex probability amplitudes c nm (t) corresponding to the Fock state with n 
photons in mode a and m photons in mode b. Assuming that the initial state is the vacuum state in both modes 
\^(t = 0)) = |0) Q |0)h these equation can be obtained in the form (we use units of h = 1): 

i-irCnm = i {n(n - \)\a + m(m - l)xt} c nm 




(4) 



Inspection of the form of Eqs. @ leads to the conclusion that the dynamics of the system is not restricted to a closed 
set of two-mode states. The number of states involved considerably in the dynamics depends straightforwardly on the 
level of excitation g - the larger the value of g the larger the number of involved states. Assuming the initial vacuum 
state (\^(t = 0)) = |0) a |0)b) and lossless dynamics only two-mode states having equal photon numbers in modes a 
and b are involved in the dynamics. 

To be more specific, the number of states that must be considered depends on the value of ratio r = g/(\ n . + Xb) 
and length of the time-evolution. For small values of r < 0.1 the conditions resemble those discussed in [40|-|42| where 
nonlinear quantum scissors have been studied. In this case, only states with small photon numbers in both modes a 
and b are involved. In particular, we can restrict the dynamics of the system to just first three lowest states of the 
whole set of two-mode states: |0) a |0)b, |l) a |l)b, |2) a |2)i>. Validity of this approximation can be judged comparing the 
approximative solution with the full (numerical) solution. Fidelity between the two states has been found useful in 
this case. A typical example investigated in Fig. 1 shows that the dynamics of the whole system can be restricted 
to the subspace of three states with accuracy of order 3 x 10~ 4 . For this case we assume that g = 0.15. This value 
seems to be the optimal one, since it is small enough to keep the system's evolution closed within the desired range 
of the states and to derive the analytical formulas for the dynamics. Simultaneously, the value g = 0.15 is sufficient 
to observe the Bell-states generation. The reason is that if we assume the value for g to be too small, mostly the 
vacuum state is populated and an efficient Bell-state generation is not possible. 

In this approximation, the evolution of probability amplitudes is governed via the equations: 

^c 00 (£) = g*cu(t) 

i^-cu(t) = gc 00 (t) + 2g*c 22 (t) 
dt 

i^-c 22 {t) = 2g Cll (t) + Ac 22 (t), (5) 
dt 

where A = Xa + Xb- 

The solution of Eqs. §5§ can be obtained analytically: 

coo(t) = r i exp (s-it) + r 02 exp (s 2 t) + r 03 exp (s 3 t) 
cii(i) = rn exp (s-it) + r 12 exp (s 2 t) + r 13 exp (s 3 t) 
c 22 {t) = r 21 exp (sit) + r 22 exp (s 2 t) + r 23 exp (s 3 t), 

(6) 
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where 
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Amplitudes take the following form: 
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Denominators occurring in Eqs. 



for amplitudes are given as: 

mi = |Ai| 2 , 

?n 2 = mt = (— 1 — 15a; 2 



M ) Ai. 



Also the following parameters have been used in the above relations: 

Ai = (3i -VE)(1+ 15a; 2 ) 4- (3? + V3)M 2 , 
X 2 = (i- V3)(l 4- 15a; 2 ) + (i + V3)M 2 , 



M = (-1 - 9x 
A 



3ia;A 



^1/3 



^3 4- 66a; 2 + 375a; 4 , 

a; = g/ixa + Xb)- 



(7) 



(8) 



(9) 



(10) 



The analytical solution in Eqs. ((BJ) valid for three two-mode states allows to analyze entanglement in the system. 
Here, we pay attention to qubit-qubit entanglement and quantify it using negativity. This entanglement measure is 
defined as 15a, k33 : 



N(p) = max ^0, —2 min/ij^ , 



(11) 



where fij are the eigenvalues of the partial transpose of the density matrix p T . The factor 2 appearing in this 
definition is chosen to get N(p) = 1 for Bell's states. Negativity is a commonly used entanglement measure allowing 
the entanglement quantification - for maximally entangled state we have N(p) = 1, whereas N(p) = occurs for the 
cases when the entanglement disappears completely. Moreover, it can be applied for more general models exceeding 
the simple qubit-qubit one. 
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In our system, we can define three subsystems of the qubit-qubit form: the first subsystem is spanned by the states 
|0) a |0)b and |l) a |l) b (negativity N i W ), the second one by states |0) a |0)b and |2) a |2) b (negativity A0220) an d the 
third one by states |l) a |l)b and |2) a |2)b (negativity iVi22i)- We note that these subsystems are not independent - 
each of the two-mode states is involved in two considered qubit-qubit subsystems. This also means that there occur 
correlations in the presence of entanglement in these subsystems. 

This is illustrated in Fig. 2, where the negativities A*oiio, A0220, and -ZV1221 for the subsystems are plotted. Negativity 
A0220 oscillates with a relatively small amplitude (around the value of 0.17) similarly as negativity -/V1221 (having the 
amplitude around 0.3). On the other hand, values of negativity iVono show that the system can be treated as a 
generator of Bell states in certain time instants in which iVoiio ~ 1- The generated Bell states are of the form: 

\B 1 ) = -L(\0) a \0) b + i\l) a \l) b ), 

\B 2 ) = -L(\0) a \0) b -i\l) a \l) b ). (12) 



Moreover, whenever negativity iVoiio achieves the shallower minimum (with amplitude approx.^ 0.3), entanglement 
is partially transferred to the subsystem spanned by states |l) a |l)& and |2) a |2)f, {1221}. In addition, when negativity 
-^0110 equals zero negativity -/V1221 also vanishes. However, the whole system is not in a pure state in this specific 
case because A0220 7^ 0. The values of negativity Aq 2 2i are not high for a fixed value of g, but we have observed that 
the mechanism of entanglement transfer alters as g increases. At larger values of g temporal evolution of negativities 
AWo j -^0220* an< ^ -^ 1221 ^ s more complex (see Fig. 3a). In particular, the role of higher populated two- mode states 
(|2) a |2)b) is more pronounced and entanglement in subsystem composed of states |0} a |0)& and |l) a |l)b is weaker. There 
occur time instants in which negativity iVono achieves one of its minima and nearly simultaneously entanglements 
in other both subspaces reach their maximal values. For comparison, we have drawn in Fig. 3 also parameter R 
that quantifies the violation of Cauchy-Schwartz inequality (CSI, for definition, see later) and serves as a measure of 
non-classicality of the whole system. 

Inspection of curves in Fig. 3 indicates that maximum of parameter R can be reached in time instants in which 
qubit-qubit negativities are far from their maximal values. This reflex complexity of quantum states found in temporal 
evolution of the discussed system. It should be noted that, for the case discussed here, we have assumed the greater 
value of parameter g. The used value g = 0.6 has been found optimal for the entanglement generation within the 
considered qubit-qubit subsystem. However, the states involving more than two photons are non-negligibly populated 
and so we cannot apply the analytical formulas derived above. 



III. ENTANGLEMENT GENERATION IN DAMPED RESERVOIR 



As any real physical system is not completely isolated, we should include into considerations about entanglement 
formation and its time-evolution also environmental effects. We describe damping in the standard Born and Markov 
approximations. We thus apply the following master equation for the reduced statistical operator p: 



dt P 



-i(Hp - pH) + la 



apa) — — [pa) a + a^ap) 



+ n<ila [o pa + apa) — a'ap — pad)] 



lb 



10 - - (0b + b r bp 



+ nblb b^ pb + b0 — b'bp — pbw , 



(13) 



where are damping constants and n aj b denote mean photon numbers of noise in modes a and b. 

Operator equation written in Eq. (|13|) can be transformed into the following differential equations of motions for 
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the elements p nm .ki of statistical operator; p nm ,ki = (k\(n\p\m)\l): 



d_ 1 

dt p nm .ki -- 2 



«Xa - 1) - m(m - 1)) 

+ i Xb (k(k - 1) - 1(1 - 1)) 

+ 7a (n + m — 2n a (n + m + 1)) 



+ 7 fc (fc + Z-2n (fc + Z + l)) 

+ gVnfc Pn-lrn,k-ll 



Pn m.k I 



-gy/(m + l)(l + l) 

+ ffV(ra + l)(fc+ 1) 
— g*\[ml pnm-l,kl-l 

+ n a y/nm p n -im-l,kl 



(1 + n b )y/ (k + + 1) p nm ,k+U+i 



■ 76 



+ n b Vkl Pnm,k-ll-1 



(14) 



In numerical calculations, the number of equations in Eq. (|14[) that have to be considered depends strongly on the 
value of constant g: the larger the value of g the greater the number of equations. 



A. Zero-temperature reservoir 



The zero-temperature reservoir is defined by the requirement that mean number of photons is zero, i.e. the system 
is influenced only by vacuum fluctuations of the environment. In the analysis we restrict ourselves to the subspace 
that contains states with no more than two photons in mode a and b. We have found for the case discussed here that 
some amount of entanglement described by negativities occurs in the system regardless on the damping strengths 
7 a . b- It is documented in Fig. 4 where the map of negativity iVi22i as a function of time t and damping constant 7 is 
plotted (7 = 7 a = 7 6 ). 

A typical temporal evolution of negativity is such that after oscillations at the beginning an asymptotic nonzero 
value is reached (set the inset in Fig. 4). It is interesting to note that the asymptotic value does not depend on 
the level of damping (for g = 0.6, A1221 ~ 0.054). On the other hand, negativity Aono shows that the effect of 
sudden death of entanglement occurs in the subsystem with zero and one photons in modes a and b. Even more 
interestingly, negativity A0220 of the subsystem composed of states |0) a , \2) a , |0) o , and |2) brings evidence of the 
effects of sudden death and sudden birth of entanglement. After entanglement rebirth negativity -ZV0220 reaches its 
nonzero asymptotic value (A0220 ~ 0.035 for g — 0.6) regardless of the damping strength 7. Positions of the borders 
between zero and non-zero values of negativities A^oiio and A0220 are plotted in Fig. 5. It is worth mentioning that 
entanglement described by negativity -/V0220 lives almost two times longer before its sudden death compared to that 
described by negativity iVono- 

Analytical treatment enables to reveal the origin of the observed features to certain extent. If the system is not 
damped, the generated states are in a coherent superposition of states |«) a |*)b containing i photon pairs [compare 
the form of Hamiltonian H lrA in Eq. |T])]. In a damped system, depopulation channels due to spontaneous emission 
are opened and, as a consequence, the system has to be described by statistical operator p. Nevertheless, projected 
statistical operators p red of the considered qubit-qubit subsystems attain an "X" type form: 



/ On 








a u 





022 














033 





V 041 








044 



(15) 



In Eq. (|15p elements an describe populations of the corresponding states whereas (i 7^ j) stand for coherences 
(probability flows). 
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The statistical operator p led written in Eq. ([15]) can be after partial transposition easily diagonalized. Its eigenvalues 
can be written as: 



Ai — an, 

a-22 + »33 - \J {a-22 - 033) 2 + 4ai 4 a 4 i 



a-22 + a.33 + V («22 - a 33 ) 2 + 4ai 4 a 4 i 



A 4 = a 44 . (16) 

Nonzero negativity N occurs whenever at least one of the eigenvalues written in Eq. (fTB|) is negative [5|| . The analysis 
of the formulas in Eq. (|16p shows that this situation is found whenever 022033 < 014041. The equality 022033 = ai 4 a 4 i 
then represents the border condition that localizes the effects of sudden death and sudden birth of entanglement. A 
detailed discussion of the relation between the form of a statistical operator and conditions for entanglement vanishing 
or rebuilding can be found in the review article Q. 

In detail, negativity iVono is larger than zero provided that 



P00,llPll,00 
P01,0lPl0.10 

Similarly, the conditon 

POO, 22/022,00 
A)2, 02/^20, 20 



< 1- (17) 



< 1 (18) 



guarantees nonzero values of negativity Nq220- I n general and roughly speaking, the conditions written in Eqs. (|17|) 
and (|18|) say that entanglement is lost at the moment when populations become larger than coherences. However, 
if coherences become larger than populations later, entanglement is revealed. This is documented in Fig. 6 for the 
evolution of entanglement in the subsystems described by negativities iVono and Nq220- Whereas the curves in Fig. 6 
clearly identify the instant of sudden death of entanglement in case of negativity iVono (Fig. 6a), instants of both the 
sudden death and sudden birth of entanglement can be determined in case of negativity ./V0220 (Fig- 6b). 

We note that the explanation of occurrence of entanglement sudden death and sudden birth presented here is very 
similar to that discussed by Ficek and Tanas in (30| when investigating two two-level atoms in a cavity. They have 
shown that dark periods of entanglement are observed provided that the population of a symmetric state became 
larger than the appropriate coherences. 

We would like to stress that, in our model, entanglement does not vanish in the long-time limit even for a damped 
system. Some parts of the whole system arc disentangled as evidenced, e.g., by zero asymptotic values of negativity 
iVono- However, there exist parts of the whole system that remain asymptotically entangled. Moreover, the study 
of temporal dynamics of entanglement has revealed that entanglement can "flow" among different subsystems. This 
may preserve entanglement from damping processes and its loss. We also note that the studied system does not evolve 
into a pure state which distinguishes it from many other systems that exhibit total disentanglement in the long-time 
limit. 



B. Thermal reservoir damping 

Considering a thermal reservoir at finite temperature, mean numbers n a and hb of reservoir photons are nonzero and 
so dynamics of the system is governed by the master equation in its general form given by Eq. (|13[) . Unfortunately, 
such reservoirs lead to weakening of entanglement compared to the case of zero-temperature reservoirs discussed 
above. 

Assuming smaller values of mean reservoir photon numbers n a and Tib, entanglement is generated in all three 
considered qubit-qubit subsystems, as described by negativities iVono, -^0220, and -/V1221 for finite times. Whereas 
entanglement indicated by negativities AToiio and -/V0220 arrives at the crucial points of sudden death, certain amount 
of asymptotic entanglement is preserved as evidenced by nonzero asymptotic values of negativity ./V1221 ■ It holds that 
the larger the mean reservoir photon- numbers n a and fib, the smaller the asymptotic values of negativity -/V1221 (see 
Fig. 7). 

The dynamics of entanglement is in general more complex for finite-temperature reservoirs compared to those at 
zero temperature, as documented in Fig. 7 where many periods of entanglement sudden death and sudden birth 
occur in the evolution of negativity iVi22i- We have also observed, that entanglement described by negativity iVono 
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disappears approximately two times faster compared to that monitored by negativity A0220 ■ The instants of sudden 
death of entanglement depend strongly on damping constants j a and 7;,: the larger the values of constants j a and 
76, the sooner the effect of sudden death occurs. Qualitatively the same dependence as that depicted in Fig. 5 for 
zero-temperature reservoirs has been revealed. 

Asymptotic values of negativity -/V1221 are zero for sufficiently large values of mean reservoir photon numbers n a 
and fib- In this case, entanglement in the whole system is completely lost. In this case, also parameter R describing 
violation of Cauchy- Schwartz inequality is smaller than one indicating classical behaviour of two modes. 

We note that, similarly as in the case of zero-temperature reservoirs, entanglement is present whenever populations 
are larger than coherences. 

Comparison with the model of a nonlinear coupler interacting with a thermal reservoir and excited by an external 
coherent field as analysed in [42[ underlines a distinguished property of this system - the ability to generate asymptotic 
entanglement at finite temperatures. 



IV. NONCLASSICAL CORRELATIONS OF INTEGRATED INTENSITIES 



Entanglement between fields in modes a and b has been monitored using negativity of different qubit-qubit sub- 
systems defined inside modes a and b. The presence of entanglement in at least one of the subsystems reflected 
entanglement between optical fields in these modes. We have observed that whenever the system shows entanglement 
also nonclassical correlations between integrated intensities of fields in modes a and b occur [5^] . These nonclassical 
correlations of intensities occur only provided that the joint quasi-distribution (in the form of a generalized function) 
of integrated intensities related to normal ordering of field operators attains negative values for some regions of inten- 
sities [4(j. These negative values can be monitored, e.g., using second-order intensity correlation functions I^ 2 ) that 
violate the Cauchy-Schwartz inequality. Violation of this inequality can be quantified using parameter R defined a 

R(t) = ■ (19) 

Vr$(*)r$(*) 

where the second-order intensity correlation functions r*- 2 -* are defined as 



r (2 
1 k 



l )(t) = Tr{p(t)Mi k I l }. (20) 



In Eq. (|20[) . p means statistical operator, Ij. = k'k and AT stands for normal-ordering operator. Nonclassical states 
obey the inequality R > 1. 

Considering the system without damping, fields in modes a and b are always entangled and also parameter R 
is always larger than one (compare temporal evolution of negativities N and parameter R in Fig. 3a). Consid- 
ering zero-temperature reservoirs, nonclassical correlations of intensities are always observed. On the other hand, 
finite-temperature reservoirs gradually destroy nonclassical correlations of intensities, similarly as they weaken en- 
tanglement. Nonclassical correlations in intensities are asymptotically preserved for smaller values of mean reservoir 
photon numbers n a and However, larger values of mean photon numbers n a and n& result in the loss of nonclassical 
correlations in finite times. We also note that the greater the values of damping constants j a and jb the sooner the 
nonclassical correlations are lost. These features are documented in Fig. 8, where temporal evolution of parameter R 
is plotted. 

Investigation of entanglement based on negativities of qubit-qubit subsystems composed of Fock states with low 
numbers need not guarantee a correct determination of non-classicality of intensity correlations. The problem is that 
entanglement "flows" among different subsystems during its temporal evolution and it may happen that it exists only 
in subspaces composed of Fock states with larger numbers. As an example we consider systems described in Fig. 8 
(n a — rib — 0.4, g = 0.6) and exhibiting entanglement also in subsystems involving states |3) a |3)f, and |4) |4)b. 



V. CONCLUSIONS 



We have analysed entanglement between number states generated in two nonlinear Kerr oscillators pumped by 
optical parametric process using qubit-qubit negativities. We have shown that maximally-entangled states of Bell 
type can be generated under certain conditions. Interaction of the system with zero-temperature reservoirs on one side 
weakens the ability to generate entanglement, but on the other side leads to a more complex evolution of entanglement 
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with the effects of its sudden death and sudden birth. The effect of sudden birth occurs provided that coherences in 
the system dominate populations. Instants of sudden deaths and sudden births of entanglement differ for different 
qubit-qubit systems which reflects "dynamical flow of entanglement" in the system. Finite reservoir temperatures 
inhibit the effect of entanglement sudden birth and destroy asymptotically entanglement for greater reservoir mean 
photon numbers. We have found that entanglement occurs whenever there exist nonclassical correlations in intensities 
of two oscillator modes that violate Cauchy-Schwartz inequality. 
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time 

FIG. 1. Declination of fidelity F from unity as a function of scaled time t. Fidelity F is calculated between the state denned in 
the subspace spanned by states |0) a |0)t,, |l) a |l)f,, and |2) a |2)f, and the state evolving in the subspace containing states |i)a|j)b 
for i,j = 0, .. ., 10. The scaled time t is measured in l/xa,b units. Initial vacuum state |0) a |0}t is assumed, g = 0.15. 
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FIG. 2. Negativities iVoiio, -^0220, and -/V1221 as functions of scaled time t measured in l/xa,b units. The coupling g = 0.15, 
other parameters and used units are the same as in caption to Fig. 1. 
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FIG. 3. (a) Negativities iVono, -^0220, and -/V1221 and (b) parameter R indicating violation of the Cauchy-Schwartz inequality 
violation as they depend of scaled time t. The scaled time t is measured in l/xa,b units. Initial vacuum state |0) a |0)(, is 
assumed, g — 0.6. 
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FIG. 4. Map of negativity N1221 as a function of scaled in l/xa,b units time t and damping 7 (measured in Xa,b units). Nonzero 
values of negativity -/V1221 occur inside grey and black areas. In the inset temporal evolution of negativity N1221 for 7 = 0.0025 
is shown. The coupling g = 0.6, other parameters and units are the same as in caption to Fig. 3. 
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FIG. 5. Borders between the regions with zero/nonzero negativities JVoiio (dashed curve) and N0220 (solid curves) at the time 
axis as they depend on damping 7 (measured in \a,b units). The curves indicate instants in which the effects of sudden death 
and sudden birth of entanglement occur. In the area between two solid lines, no entanglement described by negativity -/V0220 is 
found. The coupling g = 0.6, other parameters and the units used here are the same as for Fig. 3. 
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FIG. 6. Time evolution of quantities F (F = poo.iipn.oo) and G (G — poi.oipio.io) giving numerators and denominators, 
respectively, in conditions written in Eqs. (|17p (a) and Eqs. (|18p . Crossing of the curves indicate instants of sudden death and 
sudden birth of entanglement. Inset in (b) shows a detail of temporal evolution; 7 = 0.01. The coupling g — 0.6, other other 
parameters are the same as in caption to Fig. 3. Time t is measured in 1/Xa,b units. 
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FIG. 8. Temporal evolution of parameter R for damping constants 7 = 0.001, 0.005, and 0.01 assuming n = 0.4. For comparison, 
temporal evolution of parameter R for zero-temperature reservoirs (n = 0) is given for 7 = 0.01. The coupling g = 0.6, Values 
of other parameters are the same as in caption to Fig. 3. Time t is measured in l/xa,b units. 



